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INTRODUCTION 

The study of the anisotropies in the arrival directions of cosmic rays (CR) provides a 
handle, together with the spectrum and composition studies, to face several open ques- 
tions in the field like: which are the sources of cosmic rays?, how do CRs propagate?, 
which is the CR composition?, how are the galactic and extragalactic magnetic fields? 

The cosmic ray flux has proven to be very close to isotropic and thus very careful 
studies are needed to measure the anisotropies. Before entering the discussion of the 
techniques used, which are the main scope of these lectures, we will review some 
topics that are relevant to understand the kind of anisotropies that can be expected in 
different energy ranges. These are: the matter distribution in our neighbourhood, the 
galactic and extragalactic magnetic field effects on the CR propagation (deflections, flux 
amplification, multiple images) and the Greisen Zatsepin Kuzmin (GZK) horizon. 



Our local universe 

The distribution of matter around us is not isotropic. The Milky Way is a spiral galaxy. 
The disk, that has a width of about 1 kpc (1 pc ~ 3 ly), contains in the spiral arms 
most of the visible stars, as well as atomic gas with a density of approximately one 
particle per cm-^. The solar system is located at 8.5 kpc from the galactic center. A 
regular magnetic field component permeates the disk. The Galaxy is also permeated 
by a turbulent magnetic field where charged cosmic ray particles are trapped and that 
probably extends few kpc outside the disk. The Galaxy has also a spheroid (or stellar 
halo) of more spherical shape constituted by an older stellar population and a smaller 
amount of gas (with density ~ 0.01/cm^), that extends up to a distance of about 15 kpc. 
It has a central bulge with a bar-like shape. A larger dark matter halo is also present, 
which existence is evidenced by its gravitational effects, such as the rotation curves of 
stars and gas that remain flat up to distances larger than those containing the visible 
matter, but its composition is not yet known. The Galaxy has also a massive black hole 
in the center with a mass of few lO^M© . 



The Milky Way and its neighbour galaxy Andromeda are the most prominent mem- 
bers of a small cluster of about 30 galaxies called the Local Group, with a diameter of 
about 2 Mpc. 

The Local Group forms part of the Local Supercluster, a large agglomerate of about 
100 galaxy clusters and groups, that is dominated by the Virgo cluster located near its 
center (distant about 18 Mpc from the Milky Way). Virgo is a large cluster, with more 
than 2000 galaxies, including the prominent radio galaxy M87. 

At even larger scales the distribution of clusters shows an inhomogeneous distribution 
with large filaments and voids. Only at scales larger than few hundred Mpc the universe 
starts to look homogeneous. 



Magnetic field effects on the CR propagation 

The galactic magnetic field has a regular and a turbulent component. Despite a con- 
siderable observational effort, both are still poorly known [[l], From the observation 
of polarized radio emission of other face-on spiral galaxies it is known that the regular 
magnetic field follows the spiral pattern of the stars distribution, while radio polarization 
measurements in edge-on galaxies show the existence of magnetic halos extending few 
kpc above and below the galactic disks. 

In our own Galaxy, Faraday rotation measures of pulsars and extragalactic radio 
sources indicate that the regular magnetic field follows the spiral pattern in the disk, 
with a local value Breg — 2 — 3 /iG. [1]. The turbulent field is thought to be larger in rms 
amplitude than the regular one, with a coherence length of approximately 100 pc. 

The propagation of charged CRs is affected by the magnetic fields (both galactic and 
extragalactic) present along their trajectories. The deflection of ultra-relativistic charged 
particles is proportional to their charges and inversely proportional to their energy 

F = myv = X B, v = cu. 

Then, the direction of propagation is given by 

u = uo + ^ fdluxB. 

A CR of charge Ze and energy E in a constant field B describes a circle of Larmor 
radius 

E/Z 

Then, for E/Z < 10^^ eV (= 1 EeV) cosmic ray trajectories wind around the magnetic 
field lines and they remain confined by the galactic magnetic field for a long time, 
while for larger energies their trajectories are only moderately perturbed and they are 
not confined in the Galaxy. For E/Z < 10^^ eV they scatter off the turbulent magnetic 
field irregularities, that have a coherence scale / < r^, make a random walk and diffuse. 

Liouville theorem. An isotropic flux of CRs remains isotropic after propagating 
through a magnetic field. According to the Liouville theorem the phase space distri- 
bution f{r,p) is constant along cosmic ray trajectories if there are no processes that dis 
troy or create particles. The intensity of CRs, defined as the number of particles cross- 
ing a unit area dA per unit time dt and unit solid angle dQ. with energy within E and 
E + dE is given by / = dN/ {dAdtdQ.dE) and we can write dN = f{j,p)d^rd^p, with 



(Pr = dAvdt and d^p = p^dpdO.. Then / = f{r,p)vp^dp/dE = f{j-,p)p^. Since p is 
constant along the trajectory, / is also constant. As a consequence, an isotropic CR flux 
remains isotropic unless there is a 'shadowing effect', i.e. there exist directions from 
which particles cannot reach the detector coming from infinity. For example at low en- 
ergies this happens because of the Earth 'shadow': trajectories of antiparticles leaving 
from the detector hit the Earth due to the deflections in the geomagnetic field. This in 
particular gives rise to the so-called East- West effect: Protons with energy smaller than 
few GeV are not able to reach the Earth from the east. The sign of this E-W asymmetry 
was used to infer that CR primaries are positively charged. Towards the poles the thresh- 
old is smaller and this has the effect of increasing the CR intensity with latitude at low 
energies ('latitude effect'). 

Deflection of charged particles in the Galactic magnetic fleld. As we have discussed 
the deflection of charged particles is inversely proportional to their energy, and only for 
E/Z ^ lO^^eV the deflections in the Galactic magnetic field are expected to become 
smaller than a few degrees and CR astronomy could become feasible. If the Galactic 
B field (and composition) were known, one could correct the arrival direction to search 
for the source. This could be done by 'backtracking antiparticles' leaving the detector 
with the reverse velocity of the incoming CR. According to the Liouville theorem the 
magnetic fields cannot produce anisotropics in an isotropic flux, but they do in fact 
affect anisotropic fluxes in different ways: the flux can be amplified in some regions 
and deamplified in others in an energy dependent way, and even multiple images of 
sources located in some regions can appear. These effects can be visualized plotting for 
a regular grid of arrival directions at Earth the corresponding directions from which 
the particles arrived to the galactic halo. This is shown in Figure [T] for a particular 
model of the regular magnetic field and for particles with E/Z = 20 EeV. One may 
picture this distorted image of the sky seen from the Earth as a sheet (the 'sky sheet') 
that can be stretched and folded. A source located in a fold of this sky sheet will have 
multiple images, i.e. cosmic rays of the same energy can arrive to the Earth from several 
different directions. Moreover, the flux coming from a source in a region where the sheet 
is stretched will appear demagnified while that from a source in a compressed region will 
appear magnified 

If several CRs with different energies coming from one source are detected it would 
be possible to measure the integrated perpendicular component of the magnetic field 
along the CR trajectory and locate the actual source position. 

Magnetic lensing phenomena also appear for turbulent fields [4]. The deflections in 
a turbulent field can be viewed as a random walk of the CR particles that suffer in 
each field domain a deflection in a random direction. The accumulated deflection after 
traversing a distance L in a field with coherence length Lc is given by 

8rms = a/ L/ Lc{ZeBrmsLc/ E), 

where the last factor corresponds to the mean deflection in one domain. 

For Brms — Breg, dcflcctions in the turbulent field are smaller than those produced by 
the regular field, but can dominate the magnetic lensing effects. Multiple images appear 
below a critical energy Ec, such that typical transverse displacements among different 
paths become of order the correlation length of the B field {5rms — Lc/L). Typically 
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Figure 1. 'Sky sheet' in galactic coordinates. 

£c ^ 4 X IQi^eV Z(5„„,/5/iG)(L/2kpc)3/2(L,./50pc)-i/2. For E < , the number of 
images grows exponentially. A regime is reached with a large number of images, spread 
over a region of size and with mean magnification (A) ~ 1 (like twinkling stars). 

Magnetic fields are also present outside galaxies, but the observational constraints 
are still very poor. The amplitude in the central region of clusters may reach the /iG. 
The distribution is believed to follow the filamentary pattern of the large scale matter 
distribution. In most of the space it is usually assumed that B^ms = ~ G and 
the coherence length Lc ~Mpc. 



The Greisen-Zatsepin-Kuzmin horizon 

Soon after the discovery of the cosmic microwave background (CMB) radiation it 
was realized by Greisen, Zatsepin and Kuzmin (3] that the fluxes of CR protons with 
energies of order 10^*^ eV and above would be strongly attenuated over distances of a 
few tens of Mpc. This is due to the energy losses caused by the photo-pion production 
processes in the interactions of the protons with the CMB photons. Similarly, if CR 
sources accelerate heavy nuclei, these can photo-disintegrate into lighter ones as they 
interact with CMB and infrared (IR) photons on their journey to us. In this way the 
fragments may arrive at the Earth with significantly smaller energies than the parent 
nuclei produced at the sources. Moreover, both protons and heavy nuclei can further 
loose energy by pair production processes, although due to the small inelasticities 
involved the typical attenuation length associated with e+e^ production at ultra-high 
energies is large (~ 1 Gpc for protons). These processes limit the distance from which 
ultra-high energy cosmic rays can arrive to the Earth. For example, at energies above 
~ 60 EeV, cosmic rays should mostly come from nearby sources: in the hypothesis of an 
homogeneous distribution of sources, 90% of the flux should come from sources closer 
than ~ 200 Mpc and 50% from sources closer than ~ 100 Mpc. The distance limits turn 
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Figure 2. Exposure as a function of declination 

out to be similar for iron nuclei (although the responsible processes are different), while 
lighter nuclei should come from a much closer neighbourhood |@]. 



Exposure 

For any CR anisotropy analysis it is very important to have a good estimate of the 
expected flux in any direction for the given experimental setup in the case the CR flux is 
isotropic. The exposure measures the time integrated effective collecting area in units of 
km^ yr. For each direction of the sky (y(5, a) gives the relative exposure. For a detector 
in continuous operation it is uniform in right ascension a, being only a function of the 
declination 5. 

If the detector is fully efficient for particles arriving with zenith angle 9 < 0,„, the 
exposure has only a cos(0) modulation due to the change in the effective detection area. 
The zenith of a detector located at latitude do corresponds to a position in the celestial 
sphere given = (cos5ocosa^,cos^sina^,sin5o), where is the right ascension 
of the detector zenith. A source in a direction s— (cos 5^ cos a^, cos 5^ sin a^, sin 5^-) is 
seen at a zenithal angle cos0 = ^ ■ s = cos 5o cos 5.5 cos (of^ — a^) + sin5osin5i. The 
exposure towards a direction s is proportional to the integral of cos6{t) for all the 
times when 0{t) < 6,„, and is given by (o{5)oc cos 5 cos ^osina^ + a^ sin 5 sin where 
a,n = arccos(z) (if — 1 < z < 1) with z = (cos a,„ — sin 5o sin 5) / cos 5 cos 5o, oCm = n (if 
z < — 1) and «,„ = (if z > 1) [7]. The exposure is shown in Figure [2] for a detector at 
latitude —35°, corresponding to the Pierre Auger Observatory. 

Perfect exposure holds for a continuously operating detector at the highest energies 
where every shower triggers the detector. At lower energies the more inclined showers 
are more attenuated than the vertical ones and are less effective to trigger the detector, 
and then the exposure is no longer just proportional to cos0. There are two main 
strategies to compute the exposure in these cases: a semi-analytical method that consists 
in obtaining the zenith angle distribution from the data itself and then proceed as before 
with the integration in 6{t), or the shuffling technique that consists in simulating a large 
number of fake events using the time and zenith angle distribution of real events from 
which the isotropic expectation can be obtained. 



COSMIC RAY ANISOTROPIES 



Different signals of anisotropics are expected to appear at different energies due to the 
different sources and propagation effects involved that we will briefly discuss here. 

At the highest energies the magnetic deflections are expected to be small and the 
GZK horizon limits the distance from which CR can arrive, then we might expect to 
observe only events coming from 'nearby' sources. The expected signals are then: small 
scale clustering of events coming from the same source, correlation of events with a 
population of source candidates, intermediate scale clustering reflecting the clustering 
of local sources. Lowering the energy, the deflections increase and the GZK horizon 
also increases, then the CR flux is expected to become more isotropic, but we can still 
expect some intermediate scale clustering and correlation with the sources distribution. 
Although the distribution is expected to become more isotropic when lowering the 
energy, the increase in statistics should help to detect the smaller anisotropy signal. At 
even lower energies, a large scale anisotropy signal coming from the diffusion and drifts 
of CR in the Galactic magnetic field is expected to be seen. 

If there is a neutral component of CRs a point-like signal is expected and a correlation 
with the source population at the angular resolution scale should appear. 

In addition to these intrinsic anisotropics in the CR arrival direction, there is also a 
large scale anisotropy signal expected from the motion of the detector with respect to 
the CR rest frame, the so-called Compton Getting effect. 

We will now discuss some of the techniques that are used to measure the anisotropics 
in the CR intensity /, that is defined as the number of particles per unit solid angle that 
pass per unit time through a unit of area perpendicular to the direction of observation u. 
The differential (spectral) intensity I{E) is the intensity of particles with energy in the 
interval from E io E + dE. 



Large scale anisotropies 

The first signal at large angular scales that we look for is a dipole in some direction 
j. A dipole gives rise to an intensity I{u) = Iq + hj ■ u. The amplitude is defined as 

^ ~ ijmax I mill) / {Jmax ~l~ hnin) ~ Ii/Iq- 

The most usual methods to detect a dipole are one dimensional: they only study the 
right ascension distribution (of the full data set or in a fixed declination band). The reason 
is that some experiments cannot reliably determine the dependence of the exposure in 
the declination. 

The most standard analysis technique is the Rayleigh method, that performs an 
harmonic analysis in right ascension [jsl]. For A'^ events with right ascension a,, the 

k-order harmonic has amplitude = \J ci\ + b\ and phase = dscidin{bk/ a^), with 

at, = {2/N) Zti cos{kai) and bt, = {2/N) E^j sin(fca,-). 

The significance of a given measurement, that is the probability that an amplitude 
larger or equal than the observed arises from an isotropic data set by chance, can 
be estimated by P{> r^) = exp(— A^r^/4). According to the Central Limit Theorem, 
the sum of A'^ independent random variables xi,. ..,xj\i identically distributed with any 



probability distribution function (pdf) has a pdf approaching for large A^^ a Gaussian with 
mean equal to the sum of the means and variance equal to the sum of the variances. 
Taking as the Xj variables the right ascension coordinates of the events a,, that have a 
uniform distribution in the interval [0,2n] for an isotropic distribution of CRs, we find 
that both and are Gaussian distributed with (a^) = (bk) = and o^{ak) = o^(bk) = 
2/N. To obtain the distribution of the amplitudes r^, we have to consider that for n 
Gaussian variables ji, . . . the variable z = L/(j; ~ {yi))^/^l has a distribution 
with n degrees of freedom. Then, the variable z defined as z = {N /2){a\ + b\) = 
{N/2)rl has a distribution, P{z) = exp(— z/2)/2. Changing variable to we get 

P{rj^)drk = {l/2)exp{—Nr^/4)Nri^dri^, and integrating it above a given amplitude we 
get the advertised P{> Vk). 

The Rayleigh analysis only has information on the projection of the real dipole into 
the equatorial plane. For a full sky uniform exposure experiment, r\ = Acos(5^,p), with 
ddip the dipole direction declination, and is the right ascension of the dipole direction. 

If the exposure is not uniform or there is no full sky coverage, the relation between 
ri and the original dipole components A, = A sin 5dip and A^ = A cos 5dip in the case in 
which the exposure is independent of a, is given by [@] 

ri = 

Cl +C2Aj; 

where 

ci = L^""" d5 co(5) cos 5, C2 = f!'""' d5 co(5) cos d sin5, C3 = d5 co(5) cos^ 5. 

^min ^min ^min 

Some proposals to reconstruct the three dimensional dipole can be found in refs. 

Compton Getting ejfect. If the CR flux is isotropic in a reference system S and 
the observer is moving with respect to that coordinate system with a velocity V , he 
will measure a dipolar anisotropic flux. Let /(p, r) be the distribution function of CR 
particles in the frame 5, where it is isotropic, and f'{p' ^r') that in 5', that corresponds to 
the detector frame moving with V with respect to S. Due to Lorentz invariance /(p, r) = 
f'{p' ,r'). The momentum of particles in S' is related to that in 5 by p' = 'Yv{p — {p/u)V) 
with u the velocity of the relativistic particles. For a non-relativistic motion of the 
detector we take V and jy ^ I. Then, we can write 

The intensity can be written as I{t,E,r,h) = p^f{t,r,p), then In/ = 21np — In/ 
and dlnf/dlnp = dlnl/dlnp - 2 ~ {d\nE/d\np){d\nI/d\nE) - 2 = (1 - 
/E'^) (—7) — 2 ~ —(7+ 2) for particles with spectrum / E^. Then 

I'{E') =/(l + -(7+2)cos0). 

For example for a detector moving with V = 100 km/s and 7 ~ 3 the dipole amplitude is 
A~ 1.6x 10^3 

The orbital motion of the Earth around the Sun is expected to modulate the measured 
flux of CRs due to the Compton Getting effect. The rotation velocity of the Earth around 
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Figure 3. Amplitude and phase as a function of the energy from different CR experiments iflill . 



the Sun isV = 29.8 km/s. A vertically looking detector should see a modulation of the 
intensity with the solar time/(?) = /o(l +rcos{{t — tQ)27i:/24hs)). For every detector the 
maximum appears at a solar time to = 6 hs, this can easily be understood by thinking 
about the relative direction of the zenith and the rotation of the Earth. The amplitude 
depends on the detector's latitude and can be as large as A = {V /c) (7+ 2) ~ 5 x 10^^. A 
modulation in the solar time frequency that agrees with that expected from the Compton 
Getting effect has in fact been measured at energies around 10 TeV by the EAS-TOP 
experiment [flU] and by the Tibet Air-Shower experiment ll^. The measured amplitude 
of the solar frequency modulation has also been used to estimate the spectral index, 
obtaining 7= 3.03 ±0.55 at an energy range (6 - 40) TeV [|l3l] . in good agreement with 
the direct measurement value. 



Large scale anisotropy measurements. The good agreement with the expectations 
for the solar frequency measurements, where the expectations are well known, gives 
confidence that the measurements are reliable for the sidereal frequency analysis, that 
contains the real right ascension modulation of the CR intensity, for which the expecta- 
tions are uncertain. 

A contribution to the sidereal time frequency (or right ascension) modulation is also 
expected from the (unknown) motion of the solar system with respect to the rest frame 
of the CRs. (The solar day is a bit longer than the sidereal day: 1 year = 365.24 solar 
days = 366.24 sidereal days). 

The results of the sidereal harmonic analysis by different experiments are summarized 
in Figure [3] [ 14] . The measured amplitudes in the energy range from 10^1 to 10^4 eV are 
A ~ few X 10~^. The sidereal time modulation arises from a combination of the intrinsic 
anisotropy of the CR intensity in their own frame and the contribution from the detector 
motion. If the CR plasma were at rest with respect to an inertial system attached to the 
Galactic center, the rotational velocity of the Sun around the Galaxy of ~ 220 km/s, 
would lead to a dipole amplitude A ~ few x 10 an order of magnitude larger than the 
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Figure 4. Summary of Auger 95% CL upper bounds on the amplitude of a dipolar modulation in right 
ascension and results from previous experiments. [16] 



observed values, indicating that the CR plasma corrotate with the local stars [|l5l] . 

The transport of galactic CRs in the magnetized plasma is governed by anisotropic 
diffusion, drift and convection and detailed measurements of CR anisotropics can be 
useful to explore magnetic field characteristics and the CRs transport. In particular, 
precise measurements of the anisotropics in the knee of the spectrum may be very helpful 
to understand the origin of this feature. If the knee is due to the limit in the acceleration 
power of galactic sources, a decrease in the dipole amplitude with increasing energy is 
expected as the more isotropic extragalactic component enters in the play. Instead, if 
the knee is due to the fact that as the energy increase Galactic CRs start escaping more 
easily from the Galactic magnetic field, then an increment of the amplitude is expected 
with increasing energy as CRs flow more efficiently out of the Galaxy. Experimental 
measurements of the anisotropics cannot yet settle this point. The results at the knee and 
higher energies by different experiments are summarized in Figured [16]. 

As the energy increases and the number of events diminishes it becomes more difficult 
to measure the dipole amplitude. At energies higher than 10^^ eV AGASA has found a 
Rayleigh amplitude of 4%. More recent results from the Pierre Auger Observatory have 
put upper bounds of 1 .4% above 10^^ eV and 10% above 10^^ eV as it is shown in Figure 

m 



Small and intermediate scale clustering 

Clustering at small scales can be the clue to detect repeating sources. The amount of 
clustering (or the fraction of repeaters) gives a measure of the number of sources that 
contribute to the CRs above a given energy threshold, from which we can deduce the 
local density of sources. Clustering at intermediate angular scales contains information 
on the pattern of the distribution of the local sources. 



The autocorrelation function. The autocorrelation function is a standard technique 
to analyse the distribution of points in the sky. It measures the excess (deficit) in 
the number of pairs with respect to that expected from an isotropic distribution as a 
function of the angle. For an isotropic distribution of points on the full celestial 
sphere the expected number of pairs with angle smaller than a is given by np = 
{N{N — 1) /2) (1 — cos(a)). For partial/non-uniform sky coverage the expected number 
of pairs from an isotropic flux has to be computed simulating isotropic event realizations 
following the exposure and counting the pairs as a function of the angle. The number of 
pairs separated by less than an angle a among the N events with energy larger than a 
given threshold E is 

N /-I 

= E (1) 

,=2j=l 

where a,j is the angular separation between events i and j and is the step function. 
The chance probability for any excess of pairs at a fixed angle a and energy threshold is 
found from the fraction of simulations with a larger or equal number of pairs than what 
is found in the data at the angular scale of interest. 

The results of the autocorrelation function analysis depend on the chosen values of 
a and E (with the corresponding number of events A'^ above that energy threshold). The 
fact that the deflections expected from galactic and extragalactic magnetic fields and 
the distribution of the sources are largely unknown prevents to determine these values a 
priori. The significance of an autocorrelation signal at a given angle and energy, when 
these values have not been fixed a priori, is a delicate issue. A possible solution to 
this problem was proposed by Finley and Westerhoff iITtIi . in which a scan over the 
energy threshold and the angular separation is performed. For each value of and a, 
the fraction / of simulations having an equal or larger number of pairs than the data is 
computed. The most relevant clustering signal corresponds to the values of a and A'^ that 
have the smallest value of /, referred to as /,„,„. To establish the statistical significance 
of a given excess, it is necessary to account for the fact that the angular bins, as well 
as the energy ones, are not independent. This can be done performing a large number 
of isotropic simulations with the same number of events as the data and calculate for 
each realisation the most significant deviation Z^-^. The statistical significance of the 
deviation from isotropy is the integral of the normalised /„„„ distribution above /^f"^"- 
Then the probability that such clustering arises by chance from an isotropic distribution 
can be estimated just from the fraction of simulations having fl^^-^ < f^^'^"^- 

Although the data from a number of experiments have shown a remarkably isotropic 
distribution of arrival directions, there has been a claim of small scale clustering at en- 
ergies larger than 40 EeV by the AGASA experiment [18]. The most recently published 
analysis U9] reports 8 pairs (five doublets and a triplet) with separation smaller than 
2.5° among the 59 events with energy above 40 EeV, while 1.7 were expected from an 
isotropic flux. The probability for this excess to happen by chance was estimated to be 
less than 10^^. The significance of the AGASA clustering result was, however, subject 
of debate based on the concern that the energy threshold and angular separation were 
not fixed a priori. Tinyakov and Tkachev IlCj] computed the penalisation arising from 
making a scan in the energy threshold and obtained a probability of 3 x 10^^. Finley and 
Westerhoff [17] took also into account the penalisation for a scan in the angular scale 
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Figure 5. Autocorrelation scan 



and obtained a probability of 3.5 x 10~^. The HiRes observatory has found no significant 
clustering signal at any angular scale up to 5° for any energy threshold above 10 EeV 
ill]. A hint of correlation at scales around 25° and energies above 40 EeV, combining 
data from HiRes stereo, AGASA, Yakutsk and SUGAR experiments has been pointed 
out in ref . . 

The Auger collaboration, using the events recorded by the surface detector between 
January 2004 and August 2007, has recently reported the results of the scan performed 
in energy (above 30 EeV, with 203 events detected) and angular separation (between 1° 
and 30°) [12311 . The results are shown in the left panel of Figure [51 The most significant 
excess of pairs found is for E > 51 EeV (corresponding to 27 events) and a = 11°, with 
a chance probability oi P = 1.6 x 10~^ to arise from an isotropic distribution. Above 
this energy, a broad region of low probabilities is observed, for angular scales between 
9° and 22°. The right panel of Figure [5] shows the number of observed pairs as a function 
of the angle (dots) as well as the expectations from an isotropic distribution with the 90% 
CL bars. 



Search for point-like or extended excesses of events around a given direction in the 
sky. A neutral component of cosmic rays or a powerful source of charged particles at 
high energies could lead to a point-like excess of CRs arriving from the source direction. 
At lower energies charged particles can give rise to an extended excess of events from a 
region close to the source. Then it is important to have the tools to identify and estimate 
the significance of an excess of events from some region. The events coming from a 
source will be superimposed with that from the background, so we need to make an 
estimate of the signal and its significance. 

For any given direction the first step is to measure the observed number of events 
in a window (that can be a top-hat, Gaussian, etc) around the given direction. For 
a point-like excess (as could arise for neutral primaries) the angular resolution size 
is considered. We call this number Nqn- Then we need to estimate the background. 



For this scope we can use the detector measurements in other regions of the sky, 
Nb = aNoFF with NoFF the events measured in the OFF region and a = toN/toFF 
= (OoN / (OoFF- From this we can estimate the signal as Ns = Nqn — OlNqff- To estimate 
the significance, one possibility is to use the variance of the signal. As Nqn and Nqff are 
independent measurements, o^{Ns) = o^{Non) + OC^o^iNoFF) ■ Then there are different 
possibilities to estimate the variances. One possibility is to just consider two Poisson 
processes, then di{Ns) = \/ Nqn + OC^Nqff and the significance of the excess is given 
by 5i = (Nqn — OlNqff) / Ci {Ns) ■ Another possibility is to consider that to estimate the 
significance of an excess what we want is to asses the probability that it arises only 
from the background. Then we would take for the pdf of Nqn a Poisson distribution 
with mean equal to (Nb) and for Nqff a Poisson distribution with mean (Nb) / OC. 
Then the variance is given by (J-^{Ns) = {Nb){1 + Oc) and estimating {Nb) = toN{NoN + 

Nqff) / {toN + toFF) = {NoN +NoFF)a/ {'^ + cc),"^^ haved2(A^5) = ^^{Non +Noff) 
and 52 = {Nqn — ocNoff) / <^2{Ns)- An observed excess Ns can be said to be an 'S 
standard deviation detection'. If Nqn and Nqff are not too low (> 10), under the 
assumption that all the events come from the background ((A'5) = 0), the distribution of 
S approximates a Gaussian variable with zero mean and unit variance, and the Gaussian 
probability of S can be taken as the confidence level of the observational result. Tests 
with numerical simulations indicates that ^2 is a better estimator than 5i [24]. 

Another proposal to estimate the significance is to use a likelihood ratio method (Li- 
Ma [|24ll ). The likelihood ratio of the 'null hypothesis', corresponding here to no source 
{{Ns) = 0), and that all events are coming from the background, and the alternative tested 
hypothesis, corresponding to a non- vanishing source ((A^^) 7^ 0), is defined as 

^ L[data\null hyp) P{data\(Ns)=Q) 

L{data\alternative hyp) P{data\(Ns)y^Q) ' 



If the null hypothesis is true and Nqn, Nqff > 10, then v^21nX is Gaussian dis- 
tributed with = I. Writting 

P{data\{Ns)=0) = P{Non,Noff\{Ns) = 0, {Nb) = jf^ {Nqn + Nqff)) = 
Poisson{NoN\ {Non) = {Nb)) Poisson{NoFF\ {Nqff) = {Nb)/cc), 
P{data\{Ns) = Nqn - ccNqff, {Nb) = uNoff) = 

Poisson{NoN\ {Non) =Non) Poisson{NoFF\ {Nqff) =Noff), 
we can compute X as 
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The Li-Ma significance S = V— 21n A was shown to follow a Gaussian distribution 
better that 5i and ^2 Q. 

The search of excesses of flux from point-like and extended regions is common to 
both the Cosmic Ray and Gamma Ray astronomy fields. The selection of the ON and 
OFF is decided according to the observational data. For example, if a possible excess 
of cosmic rays from a predetermined direction in sky is tested, the rest of the sky can 
be used as the OFF region to determine the background isotropic expectation and the 
significance. If a blind search of excesses at a given angular scale is performed in the 
whole sky, proceeding in the same way for each particular direction, a distribution of 
significances is obtained. If no sources are present, they follow a Gaussian distribution 
with unit variance. 




Galactic center searches. A potential site for the acceleration of CRs in our Galaxy 
is the Galactic center region with its super massive black hole and high density of 
stars. Thus, it is an interesting region to look for an excess of events. The AGASA 
collaboration found a 4.5c7 excess {observed / expected = 506/413.6) in a 20° radius 
region close to the Galactic center region for the energy range lO^^-lO'^'^ eV. Being 
in the southern hemisphere, the Auger Observatory has a privileged view towards the 
galactic centre (GC), which passes at just 6° from the zenith at the site. For the same 
region where AGASA found the excess, and the same energy range. Auger data led 
to obs/exp = 2116/2159.6 [25], a result inconsistent with a large excess. Similarly, 
an excess reported by the SUGAR collaboration in a 5° region slightly displaced from 
the GC was not confirmed by Auger. A map of overdensity significances on 5° radius 
windows in the region around the GC is shown in fig. [6l together with the regions were 
the AGASA and SUGAR excesses were reported. The excesses present in this map are 
consistent with the expectations from fluctuations of an isotropic distribution. 



Searches for correlations with objects 

Given a population of candidate sources there are different proposed tests to search 
for a correlation with CR arrival directions. 

Cross-correlation function. This technique looks for an excess of CR separated 
by less than a given angle from any candidate source in the set with respect to the 
expectations from an isotropic CR distribution. The procedure is very similar to the 
autocorrelation analysis: we first count the number of pairs CR-objects as a function 
of the angle in the data. Then, we repeat the procedure for a large number of isotropic 
simulated data sets. To estimate the significance of any excess we compute the fraction 



of the simulations with larger number of pairs than those present in the data. This method 
was for example used to look for a possible correlation of CR with E > 10 EeV and BL 
Lacs with magnitude m < 18 at the experimental angular resolution scale (first found 
in HiRes data by Gorbunov et al. jl^). In that magnitude range there are 156 BL Lacs 
in the field of view of HiRes in the Veron Getty and Veron catalogue (27] and 271 
events have been reported by HiRes in that energy range. The number of observed pairs 
within 0.8° was 11, while only 3 were expected from an isotropic distribution. The 
fraction of isotropic simulations with a larger number of pairs is /(0.8°) = 4 x lO^'^. 
The penalization for searching at different angles and with different sets of objects is 
not included in that figure. A test of the signal using Auger data, with 1736 events with 
£■ > 10 EeV showed no evidence of excess of correlation with BL Lacs with the same 
magnitude limit in the same catalogue nIsA . 



Maximum likelihood ratio method. The idea in this method is, for a given model 
of the source distribution, to find out the values of the parameters that lead to a better 
agreement with the observational data. For example, if the model considered is that a 
fraction of GRs events come from a known population of sources and another fraction 
from an isotropically distributed background, we would say that from the set of A'^ 
measured events, Hs are source events and A^ — n., are background events. The probability 
that a background event arrives from any direction is proportional to the exposure 
(0{u). The probability that an event comes from one of the sources is peaked around 
the source direction s with some window Q{u,s) given by the angular resolution of 
the experiment (for charged particles a larger spread can be introduced to account for 
magnetic deflections). For M sources, the contribution of all the sources weighted by the 
exposure and eventually by a relative source intensity is added. For the simplest option 
of equally apparent bright sources we have 

j J 

The probability distribution for any event is P{u) = {ns/N)Q{u) + {N — n^) /N(o{u) . The 
likelihood for the set of events is L(n^) = nfLi^l";)- Then we search for the value 
of Hs that maximizes the ratio R{ns) = L{ns) /L{0), with L(0) the likelihood of the null 
hypothesis (n^ = 0). The significance can be estimated by computing R for a large set of 
isotropic simulations and counting the fraction with Rsim > Rdat- 

This method was used by HiRes to test the correlation of GRs with E > 10 EeV with 
BL Lacs with m < 18 at the resolution angular scale in their data [29]. For each event 
/ they used for Q{ui,Sj) a Gaussian centred at Sj with a dispersion equal to the angular 
resolution of that event. They found that XnR is maximized for = 8.0 corresponding to 
Ini? = 6.08. The fraction of simulations with higher XnR is / = 2 x 10 Thus, the results 
are similar in this case to those using the cross -correlation analysis. Some advantages of 
this method are that it can be adapted to give different weight to each candidate source, 
for example depending on the distance or known brightness in some band. It is also 
possible to consider different angular scales depending on the angular resolution of the 
event, or e.g. from the expected magnetic deflections in different directions. 



Binomial probability scan. For a given candidate source population, e. g. AGNs, 
galaxy clusters, radio galaxies, there are different parameters that will influence the 
correlation with events but that are difficult to fix a priori: the angular scale (magnetic 
deflections are not known), maximum distance to the objects (UHECR from distant 
sources will have their energy diminished by interactions with CMB through the GZK 
effect), energy threshold (only high energy events are expected to be correlated with 
local sources). The idea behind this method is to scan in the unknown parameters. For 
a given candidate source population, we can estimate the probability that an individual 
event from an isotropic flux has an arrival direction closer than some particular angular 
distance ^ from a member of the catalog, p, by computing the exposure-weighted 
fraction of the sky which is covered by windows of radius ^ centered on the selected 
objects. This will be a function of angular scale ^ and the maximum distance to the 
objects considered D„iax, p{^,Dmax)- For each energy threshold Emin, with A^^ events 
over the threshold, the number of events k correlated to the sources is calculated. Then, 
the probability P that k or more events of the total of events are correlated by chance 
with the selected objects is given by the cumulative binomial probability 

p=t[''j)p'i^-pr-^- 

The more significant correlation in the data set corresponds to the values ^, Emin and 
D,„ax that give rise to the smaller P value, P,,,,,,. The significance of a given correlation 
can be estimated performing a large number of isotropic simulations and under the same 
scan in ^, Emin and Dmax obtain the fraction / having a Pmin smaller than the data. 



Correlation of UHECR with AGNs. To test for possible correlations with extragalac- 
tic sources the Auger collaboration analyzed the arrival directions of the events above 
4 X 10^^ eV to look for coincidences with the positions of the known nearby (less than 
100 Mpc) active galactic nuclei from the Veron-Cetty and Veron catalog 'citeVC. The 
results of a scan over the angle Xjr between the events and the AGNs, the maximum AGN 
redshift considered Zmax and the threshold energy Ef/^ show a deep minimum in the prob- 
ability P of observing a similar or larger number of correlations arising from isotropic 
simulated data. This minimum is obtained for i// = 3.2°, Zmax = 0.017 (or maximum 
AGN distance of 7 1 Mpc) and Eth = 57 EeV (corresponding to the 27 highest energy 
events) [|30l. l3l|] . Only ~ 10^^ of the isotropic simulations have a deeper minimum under 
a similar scan. In particular, for these 27 events 20 are at less than 3.2° from an AGN 
closer than 71 Mpc, while only 6 were expected to be found by chance from an isotropic 
distribution of arrival directions. A correlation was first observed in the data obtained 
before the end of May 2006, with a very similar set of parameters, and fixing that set of 
parameters a priori the subsequent data up to August 2007 were studied, confirming the 
original correlation with more than 99% CL significance in the additional data set alone. 
This kind of prescribed test, using an independent data set and a priori fixed parameters, 
is the more safe strategy to prevent wrong claims with small statistics. It provides a clear 
way of assigning a significance to the observations, without relaying on penalizations 
for scanning. 



Figure 7. Map in galactic coordinates with the positions of the AGNs within 71 Mpc (stars) and the 27 
events with E > 57 EeV (circles of 3.2° radius). Shading indicates regions of equal exposure ||30|] . 



The map of the arrival directions and of the AGN positions is shown in Figure Ul 
A remarkable alignment of several events with the supergalactic plane (dashed line) is 
observed, and it is also worth noting that two events fall within 3.2° from Centaurus A, 
the closest active galaxy. A further interesting fact is that the energy maximizing the 
correlation with AGNs coincides with that maximizing the autocorrelation of the events 
themselves [23] and is also that for which the spectrum falls to half of the power law 
extrapolation from smaller energies 13211 . 



Log likelihood per event. The previous method cannot be applied when the candi- 
date source population is large, as the fraction of the sky covered, p, becomes of order 
unity. A more useful method in this case is to build a probability map for the expected 
arrival directions of events above a given threshold. The map can be constructed as fol- 
lows, a Gaussian of given o is taken around the direction of each object in the catalog, 
weighting them by a factor 

Andl{z)(l){z) JEi{z,E,h) 

where di is the distance to the object, 4>(z) is the selection function of the catalog and 
the integral term measures the fraction of the flux from a source at redshift z that reaches 
the Earth with an energy larger than the threshold Efh. Ej is the initial energy that the 
particle needs to have at the source to arrive at Earth with E^h, s is the source spectral 
index. As an example, the left panel of Figure [8] shows the map corresponding to the 
AGNs in the Veron-Cetty and Veron catalog and an energy threshold of 80 EeV. The 
likelihood associated to a given set of A'^ observed events is L = Y[f:=iP{'^i)i with P{ui) 
proportional to the map density in each event direction. In order that the mean value be 
independent of the total number of events, it is more convenient to use the log likelihood 



VC 2deg density map 




Figure 8. Likelihood map (left panel) and Log likelihood per event histogra?nstrigiit'pane15' 



per event 

i=i 

The idea is to measure LL of the data using a model reference map. Then simulate events 
distributed according to some alternative hypothesis: isotropic, following AGNs,... and 
compute LL for the reference model map. Then plot the histogram the LLs for each 
hypothesis: the mean of the distribution is independent of the number A^^ of observed 
events, but the width becomes smaller as grows. When the histograms corresponding 
to the different hypothesis do not overlap the test is good to discriminate among them. 

In the right panel of Figure [8] we show the histograms for 20 simulated isotropic 
events and for 20 events following the distribution of AGNs in the Veron-Cetty and 
Veron catalog using a Gaussian window of 2 degrees size. 

FINAL REMARKS 

A variety of different methods are needed to study anisotropics at different energies 
and angular scales. It is impossible to review in a limited space all the interesting ideas 
that have been put forward to this scope, I have only presented a (personally biased) 
subset of them, as well as a selection of the experimental results related to the techniques 
described. 

This is a very special time for the field as the CR astronomy is finally starting and 
we are getting the first clues on the UHECR origin: they are correlated with nearby 
extragalactic matter. There are still many open questions as which are the sources, which 
is the CR composition, how are the relevant magnetic fields? More data is eagerly waited 
to clarify these issues. 
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